GAS CLUSTER ION BEAM SIZE DIAGNOSTICS AND WORKPIECE 

PROCESSING 

CROSS REFERENCE TO RELATED PATENT APPLICATIONS 

1 This application claims priority of U. S. Provisional 

Application Serial No. 60/218,248 entitled CLUSTER SIZE 
DETERMINATION, filed July 14, 2000, the provisional 
application being incorporated herein by reference. 

BACKGROUND OF THE INVENTION 

*Q 2 This invention relates generally to the measurement of 

:^ the size distribution of gas clusters in gas cluster ion 

,p beams, and, more particularly to apparatus and methods for 

;W the measurement of the size distribution of gas clusters by 

*r§ making time-of -flight measurements. 

e 3 The use of a gas cluster ion beam (GCIB) for etching, 

H cleaning, and smoothing of 'material surfaces is known (see 

j«* for example, US patent 5,814,194, Deguchi et al . ) in the art. 

fy For purposes of this discussion, gas clusters are nano-sized 

1^ aggregates of materials that are gaseous under conditions of 

standard temperature and pressure. Such clusters typically 
consist of aggregates of from a few to several thousand 
molecules loosely bound to form the cluster. Such clusters 
can be ionized by electron bombardment or other means, 
permitting them to be formed into directed beams of known and 
controllable energy. The larger sized clusters are the most 
useful because of their ability to carry substantial energy 
per cluster ion, while yet having only modest energy per 
molecule. The clusters disintegrate on impact, with each 
individual molecule carrying only a small fraction of the 
total cluster energy. Consequently the impact effects of 
large clusters are substantial, but are limited to a very 
shallow surface region. This makes ion clusters effective 
for a variety of surface modification processes, without the 
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tendency to produce deeper subsurface damage characteristic 
of monomer ion beam processing. 

4 Means for creation of and acceleration of such GCIB's 
are described in the reference (US 5,814,194) previously 
cited. Presently available ion cluster sources produce 
clusters ions having a wide distribution of sizes, n (where n 
= the number of molecules in each cluster - in the case of 
monatomic gases like argon, an atom of the monatomic gas will 
be referred to as a molecule and an ionized atom of such a 
monatomic gas will be referred to as a molecular ion - or 
simply a monomer ion - throughout this discussion) . 

5 To a first order approximation, the surface modification 
effects of an energetic cluster is dependent on the energy of 
the cluster. However, second order effects are dependent on 
the velocity of the cluster which is dependent on both the 
energy of the cluster and it's mass (and hence the cluster 
size, n.) In order to maximize the utility of a gas cluster 
ion beam for surface processing, it is useful to know and 
control both the energy of the cluster and the cluster size 
or cluster size distribution. In certain applications gas 
cluster ion beams are used for deposition or growth of 
surface films. When thus used, it is important to know the 
mass flow to the workpiece. The quantity of clusters is 
readily determined by measuring the cluster ion current that 
reaches the workpiece. In the usual case, ionized clusters 
from a practical ionized cluster source, do not necessarily 
all carry the same number of electrical charges. By suitable 
selection of ionization conditions, it can be arranged that 
the cluster ions predominately carry a single electrical , 
charge, and in such case it is accurately assumed that each 
charge corresponds to a single cluster, but unless the 
average size or size distribution (average mass or mass 
distribution) is also known, the total mass flow to the 
workpiece is not known. It is also possible, by controlling 
the source conditions, to influence both the ratio of cluster 
ions to molecular ions and to influence the cluster size 
distribution as well. However, unless a means is available 
to measure and monitor the cluster size distribution and the 




ratio of cluster ions to molecular ions, meaningful 
adjustment and control of the source conditions for 
influencing ionized cluster size is difficult. When the 
ionized clusters do not all predominately carry a single 
charge or a known number of charges, knowledge of the ionized 
clusters ' mass per charge can also be used as a useful 
parameter to control or adjust the beam's effectiveness for 
smoothing, etching, or other processing. For these and 
other reasons it is useful to have a means of measurement 
that can provide cluster size distribution information about 
a gas cluster ion beam or that can provide information about 
a gas cluster ion beam's cluster-size-per-charge distribution 
or cluster-mass-per-charge distribution . 

Because molecular ions, as well as cluster ions, are 
produced by presently available cluster ion beam sources, 
those molecular ions are accelerated and transported to the 
workpiece being processed along with the cluster ions. 
Molecular ions, having high energy with low mass results in 
high velocities, which allow the light molecular ions to 
penetrate the surface and produce deep damage which is likely 
to be detrimental to the process. Such sub-surface ion 
damage is well established and well known from the more 
traditional monomer ion beam processing art and can produce a 
variety of deep damage and undesirable implantation. 

It has become known in the ion cluster beam art that 
many GCIB processes benefit from incorporating means within 
GCIB processing equipment for eliminating molecular ions from 
the ion cluster beams. Electrostatic (see US patent 
4,737,637, Knauer) and electromagnetic (see Japanese laid 
open application JP 03-245523, cited as prior art in US 
patent 5,185,2 87, Aoyagi et al . ) mass analyzers have been 
employed to remove light ions from the beam of heavier 
cluster ions. Electrostatic and electromagnetic mass 
analyzers have also been employed to select ion clusters 
having a narrow range of ion masses from a beam containing a 
wider distribution of masses. 

Presently practical GCIB sources produce a broad 
distribution of ion cluster' sizes with limited cluster ion 



currents available. Therefore it is not practical to perform 
GCIB processing by selecting a single cluster size or a 
narrow range of cluster sizes the available fluence of 
such a beam is too low for productive processing. It is 
preferred to eliminate only the molecular ions and other 
lowest mass ions from the beam and use all remaining heavier 
ions for processing. Practical experience has shown that it 
is often sufficient to provide filtering to eliminate 
molecular ions while depending on the typical cluster size 
distribution characteristics (few small size clusters are 
created by typical sources) to limit the small clusters (n = 
2 to -10) in the beam. Clusters of size n > 10 are 
adequately large to provide acceptable results in most 
processes. Since the typical cluster distribution contains 
clusters of up to n = several thousand and there are few 
clusters of mass less than 100, it is not significantly 
detrimental if clusters up to size 100 are removed from the 
beam in the process of eliminating the molecular ions. 
However, in order to adequately predict the processing 
effectiveness of a gas cluster ion beam, it is very desirable 
to know the distribution of masses or cluster sizes in the 
beam and to know whether molecular ions and the smallest size 
cluster ions are present or not. 

It is therefore an object of this invention to provide a 
way of measuring the mean cluster ion size or mass in a GCIB. 

It is also an object of the invention to provide means 
to determine the cluster size (or mass) distribution or the 
mass flow of cluster ions in a GCIB processing system without 
necessitating the rejection of a portion of the beam through 
magnetic or electrostatic mass analysis. 

It is also an object of the invention to determine the 
presence or absence in a GCIB of undesirable molecular ions. 

SUMMARY OF THE INVENTION 

The objects set forth above as well as further and other 
objects and advantages of the present invention are achieved 
by the embodiments of the invention described hereinbelow. 
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In the present invention it is preferably arranged that 
the molecular and cluster ions produced in an ionization 
system for gas clusters predominantly carry a single 
electrical charge and such ions are accelerated through a 
known acceleration potential. Based thereon, the ions, both 
molecular ions and cluster ions, in a GCIB produced in the 
apparatus of the invention have known and controllable 
energies per ion, regardless of cluster size or whether the 
ion is a cluster or a molecular ion. Since the ions are all 
generated from a pure gas, the molecular weight of each 
molecule is known. Therefore, by measuring the times-of - 
flight of the known ions of known energy over a known 
distance, it is possible according to the invention to 
calculate the cluster size distribution function for the ions 
in the beam. Alternatively, when it is not practical or not 
desired to arrange that all molecular and cluster ions 
produced in the gas cluster ionization system predominately 
carry a single electrical charge, in the present invention, 
by measuring the times-of -flight of the ions to calculate the 
cluster size per charge of the ions. The times-of -flight are 
determined by analyzing the fall-off of the GCIB current when 
the GCIB is abruptly terminated. 

By providing a suitable beam gating method and 
introducing appropriate sensing apparatus in a conventional 
GCIB processor, the invention makes an in situ capability for 
measurement of the mass and cluster size distribution of the 
ions in a GCIB. The invention can measure either the mean 
cluster size (or mass) or the distribution of cluster sizes 
(or masses) in a GCIB. 

For a better understanding of the present invention, 
together with other and further objects thereof, reference is 
made to the accompanying drawings and detailed description 
and its scope will be pointed out in the appended claims. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a graph showing typical ion cluster size 
distribution for GICB source; 
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17 Figure 2 is a schematic showing the basic elements of a 
prior art GCIB processing system; 

18 Figure 3 shows a schematic of the concept for in situ 
cluster size determination by time-of- flight measurement ; 

19 Figure 4 is a schematic showing a GCIB apparatus having 
improvements according to the invention; 

2 0 Figure 5 shows a graph of the simulated current shut off 

for a beam with N=3000 molecules and a = 400 molecules; 

21 Figure 6 shows a graph of the simulated current shut off 
for a beam with N=3 000 molecules and a =900 molecules; 

22 Figure 7 shows a graph of the simulated current shut off 
i* 5 *. for a Beam with N=1500 molecules and a =400 molecules; 

'™ 23 Figure 8 is a schematic showing an electronic beam gate 

(3 for gating a GCIB according to the invention; 

P 

24 Figure 9 is a schematic showing an extraction switching 

{<* technique for gating a GCIB according to the invention; 

eM 25 Figure 10 is a schematic of a technique for ionization 

O switching of a GCIB according to the invention; , 

~? 2 6 Figure 11 is a schematic showing a GCIB apparatus with 

y deflection switching of a GCIB according to an embodiment of 

:J the invention; 

27 Figure 12 shows a graph of a normalized negative 

derivative divided by time for Gaussian cluster distribution 
(N=1500 ,*a =400) ; 

2 8 Figure 13 shows a graph comparing the initial (Gaussian) 

distribution with derived distribution (N=1500 ; a =400 ) ; 

2 9 Figure 14 shows a graph of the current shutoff for a 

Poisson cluster distribution (N=1500; m=2 ) ; 

3 0 Figure 15 shows a graph of the normalized negative 

derivative divided by time for a Poisson cluster distribution 
(N=1500; m=2) ; and 
31 Figure 16 shows a graph comparing the initial (Poisson) 

distribution with derived distribution (N=1500 ;m=2 ) . 



Figure 17 shows a graph of the measured current shut off 
for a GCIB gated off using deflection switching according to 
the invention. 

Figure 18 shows a graph of the time derivative of the 
measured current shut off shown in Figure 17 . 

Figure 19 shows a graph of the cluster probability 
distribution function, f(n/q), corresponding to the current 
shut off characteristic measurement shown in Figure 17 and 
obtained according to the invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Figure 1 shows the typical ion cluster size distribution 
produced by a commonly available GICB source. The cluster 
formation process has been shown (in N. Kofuji, et al . , 
Development of gas cluster source and its characteristics, 
Proc . 14th Symp. on Ion Sources and Ion-Assisted Technology, 
Tokyo (1991) p. 15) to produce few small size clusters 
(values of N from 2 to about 10), but molecular ions (n = 1) 
are produced in abundance as are larger clusters ( n > a few 
tens, up to a few thousands.) It is known from the teachings 
of US patent 5,459,32 6, Yamada, that molecules in a cluster 
are not individually energetic enough (on the order of a few 
electron volts) to significantly penetrate a surface to cause 
the residual surface damage typically associated with the 
other types of ion beam processing in which individual 
molecular ions may have energies on the order of thousands of 
electron volts. Nevertheless, according to Yamada and 
Matsuo, Cluster ion beam processing, Matl . Science in 
Semiconductor Processing- I, (1998) pp. 27-41, the clusters 
themselves can be made sufficiently energetic (some thousands 
of electron volts) , to effectively etch, smooth or clean 
surfaces. 

Figure 2 shows a typical configuration for a GCIB 
processor 100 of a form known in prior art, and which may be 
described as follows: a vacuum vessel 102 is divided into 
three communicating chambers, a source chamber 104, an 
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ionization/acceleration chamber 106, and a processing chamber 
108. The three chambers are evacuated to suitable operating 
pressures by vacuum pumping systems 146a, 146b, and 146c, 
respectively. A condensable source gas 112 (for example 
argon or N 2 ) is admitted under pressure through gas feed tube 
114 to stagnation chamber 116 and is ejected into the 
substantially lower pressure vacuum through a properly shaped 
nozzle 110. A supersonic gas jet 118 results. Cooling, 
which . results from the expansion in the jet, causes a portion 
of the gas jet 118 to condense into clusters, each consisting 
of from a few to several thousand weakly bound molecules. A 
gas skimmer aperture 12 0 separates the gas products that have 
not been formed into a cluster jet from the cluster jet so as 
to minimize pressure in the downstream regions where such 
higher pressures would be detrimental (e.g., ionizer 122, 
high voltage electrodes 12 6, and process chamber 108) . 
Suitable condensable source gases 112 include, but are not 
limited to argon, nitrogen, carbon dioxide, oxygen, and other 
gases . 

After the supersonic gas jet 118 containing gas clusters 
has been formed, the clusters are ionized in an ionizer 122 . 
The ionizer 122 is typically an electron impact ionizer that 
produces thermoelectrons from one or more incandescent 
filaments 124 and accelerates and directs the electrons 
causing them to collide with the gas clusters in the gas jet 
118, where the jet passes through the ionizer 122. The 
electron impact ejects electrons from the clusters, causing a 
portion the clusters to become positively ionized. A set of 
suitably biased high voltage electrodes 126 extracts the 
cluster ions from the ionizer, forming a beam, then 
accelerates them to a desired energy (typically from 1 keV to 
several tens of keV) and focuses them to form a GCIB 12 8 
having an initial trajectory 154. Filament power supply 136 
provides voltage V F to heat the ionizer filament 124. Anode 
power supply 134 provides voltage V A to accelerate 
thermoelectrons emitted from filament 124 to cause them to 
bombard the cluster containing gas jet 118 to produce ions. 
Extraction power supply 138 provides voltage V E to bias a 



high voltage electrode to extract ions from the ionizing 
region of ionizer 122 and to form a GCIB 128. Accelerator 
power supply 140 provides voltage V Acc to bias a high voltage 
electrode with respect to the ionizer 122 so as to result in 
a total GCIB acceleration energy equal to V Acc . One or more 
lens power supplies (142 and 144 shown for example) may be 
provided to bias high voltage electrodes with potentials (V L1 
and V L2 for example) to focus the GCIB 128. 

A workpiece 152, which may be a semiconductor wafer, a 
thin film on a substrate, or other workpiece to be processed 
by GCIB processing, is held on a workpiece holder 150, 
disposed in the path of the GCIB 128. Since most applications 
contemplate the processing of large workpieces with spatially 
uniform results, a scanning system is desirable to uniformly 
scan the GCIB 12 8 across large areas to produce spatially 
homogeneous results. Two pairs of orthogonally oriented 
electrostatic scan plates 130 and 132 can be utilized to 
produce a raster or other scanning pattern across the desired 
processing area. When beam scanning is performed, the GCIB 
128 is converted into a scanned GCIB 148, which scans the 
entire surface of workpiece 152 . 

Figure 3 shows some basic concepts of the invention for 
measuring cluster size. A beam gating device (shown here as 
a mechanical shutter, for example, and different embodiments 
of which are discussed in more detail later) rapidly shuts 
off or restores the beam. When the beam is shut off by the 
beam gate, this leaves, in space, a cylinder of ionized beam, 
at an upstream location (for example at the exit of the 
ionizer) , which continues to move toward the workpiece where 
it may be measured by a current measuring device such as a 
Faraday cup dosimeter. The current to the Faraday cup would 
then be expected to remain constant until the lightest 
clusters travel the distance, L, from the shutter to Faraday 
cup and then to drop to zero by the time the heaviest 
clusters have traveled this same distance, L. The current 
measured at the Faraday cup falls to half its initial value 
in the time required for the average sized cluster to 
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transverse distance, L. In other words, the half -current 
delay time is given in Eqn. 1: 



m_ m a N 

40 T ^ L fe: 



41 where N is the average cluster size, ri\> is the mass of a 

proton, and ma is the mass in atomic mass units of a single 
molecule in the cluster, q is the ion charge (= e, the 
electronic charge) , and V Ac c is the acceleration voltage for 
the ions. Of course, Eqn. 1 can be inverted to determine N: 



2qV Acc fT 50% 



42 N = n Acc Eqn. 2 



m p m a 



43 To explore this further a simple model was developed by- 

assuming that the cluster number distribution is distributed 
according to a normal (Gaussian) distribution. In that case 
the cluster size distribution, f (n) , is 



j (n-N)" 



44 f(n) = — 2o ~ Eqn. 3 



45 where a is the standard deviation for the distribution. We 
must assure that there are no clusters with negative numbers 
of molecule in the clusters. Thus, we require 

46 N-3<7»0 Eqn. 4 

47 In terms of the current carried, although the lighter 
clusters travel faster than the heavier clusters, the lower 
density of lighter clusters in the column offsets this 
difference. Thus, in this model, the current carried by a 
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particular size cluster is proportional to the Gaussian 
probability. Current drops in steps as each cluster mass 
travels the path, L, and passes totally into the Faraday. 

According to this invention, an improved GCIB apparatus 
as shown in Figure 4 has a stationary GCIB and a mechanically 
scanned workpiece holder for scanning the wafer relative to 
the beam. Referring to the GCIB processing apparatus 2 00 of 
the invention shown in Figure 4, a substantially stationary 
GCIB 2 02 is directed at a workpiece 210 that may be a 
semiconductor wafer, a metal substrate, or the like. 
Workpiece 210 is held on an X-Y positioning table 204 that is 
operable to move the workpiece in two axes, effectively 
scanning the workpiece 210 relative to the GCIB 202. The 
GCIB 202 impacts the workpiece 210 at a projected impact 
region 230 on a surface of the workpiece 210. By X-Y motion, 
the table 204 can position each part of a surface of the 
workpiece 210 in the path of GCIB 202 so that every region of 
a surface of workpiece 210 may be made to coincide with the 
projected impact region 230 for processing by the GCIB. An 
X-Y controller 216 provides electrical signals to the X-Y 
positioning table 204 through electrical cable 218 for 
controlling the position and velocity in each of X-axis and 
Y-axis directions. The X-Y controller 216 may receive 
control signals from an overall system controller, not shown. 
X-Y positioning table 2 04 moves by continuous motion or by 
stepwise motion according to conventional X-Y table 
positioning technology to position different regions of the 
workpiece 210 within the projected impact region 230. The 
workpiece holding surface 206 of positioning table 204 is 
electrically conductive. A cylindrical enclosure 220 
surrounding the GCIB 202 upstream of the workpiece 210 is 
electrically conductive. Together the workpiece holding 
surface 206 and the cylindrical enclosure effectively 
comprise a form of Faraday cup assembly 232 for collecting 
the GCIB current arriving at the workpiece 210 and workpiece 
holding surface 206. The workpiece holding surface 206 and 
the cylindrical enclosure 220 are electrically connected 
together by electrical lead 222. An electrically insulating 
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layer 208 of positioning table 204 isolates the workpiece 210 
and workpiece holding surface 206 from the other portions of 
the positioning table 204. Electrical charge induced in the 
Faraday cup assembly by the GC1B 2 02 striking workpiece 210 
is conducted through electrical lead 222 to cluster size 
measurement system 226 for measurement. An ionized beam 
switching beam gate 212 is disposed in the path of GCIB 202. 
Beam gate 212 has an open state for allowing the passage of 
ionized beam and a closed state for stopping the flow of 
ionized beam. A control cable 224 conducts control signals 
from cluster size measurement system 226 to beam gate 212, 
said control signals controllably switching beam gate 212 to 
either of its open or closed states for enabling or disabling 
the passage, of ionized GCIB to the Faraday cup assembly 232. 
A flight distance 228 is symbolized by the symbol, L, and 
represents the distance an ion flies between the effective 
portion of the beam gate 212 and the effective location of 
the Faraday cup assembly 232. A beam filter 214 is, 
preferably but not necessarily, provided to eliminate 
molecular ions, or molecular ions and light cluster ions, 
from the GCIB 202. The beam filter uses conventional 
technology as described for example in US patent 5,185,272 or 
US patent 4,737,637. 

It is recognized that the cluster size measurement 
system 22 6 may be implemented in a variety of ways and may 
include as a part of its control systems all or part of a 
programmable computer system. The X-Y controller 216 may 
include as part of its logic all or part of a programmable 
computer system. Some or all of the X-Y controller 216 and 
cluster size measurement system 226 logic may be performed by 
a small general purpose computer (not specifically shown) 
that also controls other portions of the GCIB processing 
apparatus 200. In operation, a GCIB is established and the 
cluster size measurement system 22 6 sends signals commanding 
the abrupt closing of the beam gate 1212 to abruptly shut off 
flow of ionized beam to the Faraday cup assembly 232. The 
cluster size measurement system 226 measures the time 
variation of the GCIB current, I, collected by the Faraday 
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cup assembly 232 during time from the closing of the beam 
gate 212 until beam current has ceased to be collected by the 
Faraday cup assembly 232, this includes the flight times over 
the distance, L, of all of the (fastest through slowest) 
ions. The cluster size measurement system 226 processes the 
current and time information to measure the average cluster 
size and/or mass and/or to measure the cluster size and/or 
mass distribution in the GCIB. The processing and 
measurements are made according to other principles of the 
invention taught hereinafter. A display system 240 receives 
processed information from measurement system 22 6 and 
displays digital and/or graphical cluster size and cluster 
i~= size distribution characteristics. 

*P 50 Although Figure 4 shows a Faraday cup assembly 232 as 

i=J the means for current collection as an example and not for 

;p limitation, it is recognized that a variety of alternative 

jfcf current collection arrangements can be utilized in the 

=0 - invention. Although a beam gate 212 is shown in a particular 
! location downstream of the high voltage electrodes 126, and 

l z 

IT* is shown, for example and not for limitation, to be of a type 

I 4 similar to the repeller electrode gate shown in more detail 

:ff in Figure 8, it is recognized that a variety of means and 

1^ locations for abruptly gating off and on the flow of ions are 

allowable for effective operation of the invention (several 
are discussed hereinafter) . It is important for the 
invention to have an upstream means of abruptly gating off 
and on the flow of ions and a downstream means of measuring 
the time relationship of ion current, the gating and 
measurement separated by an ion flight distance, L. 
51 Figure 5 shows results from a computer model simulating 

the current shut off for a 20 keV beam with a distribution 
having a peak (average) of 3000 argon molecules per cluster 
and a standard deviation of 400 molecules. The 50% point 
corresponds very closely to the transit time for the peak of 
the distribution, which for the symmetrical Gaussian, is the 
average. The 10 to 90% fall time for the current is 
dependent on the width of the distribution. Figure 6 shows 
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the current shut off for a beam similar to that in Figure 5, 
but with a 900 molecule standard deviation. The slower fall 
time is evident. However, note that the distribution about 
the 50% line is not symmetric, showing that a simple time 
derivative will not return the initial distribution. 

Of course, with changes in energy and mean cluster size 
the time to the 50% point will change. Figure 7 shows the 
current shut off for a 20 keV beam with an average of 1500 
molecules per cluster and again a standard deviation of 400 
molecules. In this case the 50% point has moved in to about 
67 /isec. 

Figures 5 through 7 help to understand the requirements 
for the beam gate to be used in the invention. When we are 
primarily interested in the average cluster size then the 
shutter closure must be fast compared to the 50% delay time. 
Such a beam gate requires a switching time, for the example, 
of less than 5 ^isec. In order to extract information on the 
cluster size distribution, then the beam gate switching time 
needs to be fast compared to the current fall time. In this 
case, the closure time should be 1 Msec or less, for the 
example . 

Fig. 8 shows the design for a true electronic beam gate 
A charged beam repeller electrode is added in the beam path 
to serve as a beam gate. The resistor normally grounds this 
electrode. However, when the spark gap is triggered, the 
0.003 mF capacitor will discharge into the 10 kohm resistor 
bringing the electrode up to +3 0 kV (chosen to be greater 
than the 20 kV beam acceleration voltage, V AC c - This will 
repel the ionized beam and prevent it from going through the 
electrode opening. The capacitor will then discharge with 
the RC time constant (a few 100 microseconds) and eventually 
will recharge through the 1 Mohm resistor allowing a 
repetition . 

It is recognized that the repeller electrode shown in 
Figure 8 could be a separate electrode as shown, or could 
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even be incorporated into the set of high voltage electrodes, 
126 shown in Figure 4. In the latter arrangement, the 10 
kohm resistor would return to a power supply instead of to 
ground. Some damping in the discharge circuit can be added 
if desired to minimize ringing of the voltage on the repeller 
electrode. This damping resistance can be low since the 
stray capacitance and inductance of the circuit can be made 
quite low. In any event a 1 Msec fall time is easily 
achievable. Conventional spark gaps can easily handle this 
service (~ 1 Joule of stored energy) and can be repetitively 
pulse at a few hertz. The required spark gap is small and be 
commercially purchased at modest cost. Another concern is 
the 0.003 mF capacitor, which is can be somewhat large and 
about as expensive as the spark gap, but is commercially 
available . 

An alternative way of gating the ionized beam is to shut 
off or enable extraction of the beam from the ionizer. An 
example of this is shown in Figure 9. In this case, a 
hydrogen thyratron is used to shut off the extraction 
voltage. The stored energy is even less than in the case of 
the repeller electrode design and the capacitor requirements 
are less demanding. A conventional "doorknob" capacitor is 
readily available and can be used. A glass hydrogen 
thyratron costs a few tens of dollars and can be switched at 
10' s of Hertz with simple solid state trigger electronics. 
Stray capacitance and inductance would be more significant 
than in the repeller electrode case, but ringing prevention 
can be achieved by conventional methods. A bigger problem is 
the fact that the supplies floating at extraction voltage 
experience large transients, which could reduce their 
reliability. An alternative to extraction switching would be 
to include an additional electrode at the exit of the ionizer 
and to pulse this positively. 

Figure 10 shows an example of gating the ionized beam by 
switching the ionizer on and off. Here the anode bias 
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voltage of the ionizer, V b ias/ (corresponds to V A in Figures 2 
and 4). is used to extract electrons from the thermionic 
filaments is switched to the negative filament potential 
(negative side of Vfi lament in Figure 10, and corresponding to 
the negative side of V f in Figures 2 and 4) . This 
effectively terminates any electron emission current from 
these filaments and will shut down ionization of the clusters 
very rapidly. However, switching time is limited by the 
transit time of the ionized clusters down the (for example, 
1.5") length of the ionizer 122. This will give a switching 
time of about 5 /isec, adequate for average mass determination 
but marginal for examining the size distribution. The 
circuitry is floating at acceleration potential, so the 
insulated gate bipolar transistor (IGBT) is triggered from 
ground potential using, for example, the 500 pF blocking 
capacitor shown in Figure 10. Alternatively, a 30 kV 
isolated pulse transformer or other conventional methods can 
be used. A pulse transformer has the advantage of isolating 
the IGBT from any extraction voltage turn on transients or 
transients due to arcing in the extraction column. Note that 
this same circuitry could be used for pulsing the extra 
electrode mentioned in the previous section. A few hundred 
volts would be adequate to block all ions. 

Figure 11 shows a GCIB processing apparatus 300 
according to an alternate embodiment. In the processing 
apparatus 3 00, gating the ionized beam is accomplished by 
deflection switching the GCIB 128. In this embodiment, a 
pair of deflection plates 302 receive a high voltage 
deflection pulse 308, creating an electric field that 
deflects the gas cluster ion beam away from its initial 
trajectory. The deflection pulse 308 is generated by the 
measurement system 22 6 and has a duration, T pw , greater than 
the time-of -flight of the slowest cluster ions over the 
flight distance 228, L . The deflection pulse 308 has an 
amplitude of -Va that is typically on the order of a few 
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kilovolts and suffices to deflect the GCIB 128 from its 
initial trajectory 304 to a second trajectory 306. The 
GCIB's initial trajectory 304 strikes the workpiece 210 at a 
projected impact region 23 0 on the surface of the workpiece 
210 for processing the workpiece 210. The projected impact 
area 230 is within a Faraday cup assembly 232 for beam 
current collection. When the measurement system 22 6 actuates 
the deflector with deflection pulse 308 at a time tO, the 
initial trajectory 304 of the GCIB is deflected to the second 
trajectory 306. The second trajectory 306 impacts a beam 
stopper 308 that is outside of the Faraday cup assembly 232. 
The rise and fall times of the deflection pulse 3 08 are short 
(of the order of 1 microsecond or less) so as to abruptly 
deflect the GCIB out of or into the Faraday cup assembly 232, 
which acts as a current sensor. The current sensed in the 
Faraday cup assembly 2 32 is conducted to the measurement 
system 226 on electrical lead 222. When the deflector 
switches the beam out of or into the Faraday cup, the 
corresponding fall or rise in current collected by the 
Faraday cup assembly 232 is processed by the measurement 
system according to the method of the invention and GCIB 
cluster size or cluster size distribution information is 
displayed on the display system 240. 

There are other means for switching the beam current off 
or on quickly, but one, particularly deserves some 
consideration. Arcing normally occurs occasionally in the 
extraction electrode set. This causes a momentary drop out 
of the current in the beam. By detecting these naturally 
occurring beam dropouts and measuring the beam current shut 
off transient, the cluster size can be estimated from time- 
of -flight. The timer for the current fall would be triggered 
from the arc and the time to the drop out at the Faraday or 
other current sensor measured. 

A number of techniques can be used to shut down the 
ionized cluster ion beam and determine the effect on the 
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measured current. Gating techniques have been given for 
example, not limitation, and numerous approaches to abrupt 
gating of the beam will occur to those skilled in the arts. 
The Faraday cup or other beam sensing system, which has been 
designed for the mechanically scanned smoother must be fast 
enough to accurately record the signals shown in Figures 5-7, 
This can be achieved with conventional designs. The most 
difficult element is the beam switching, with several 
practical methods available. Ionization switching is -the 
preferred approach due to its simplicity, but the active gate 
using a repeller electrode (Figure 8) and the beam deflection 
method (Figure 11) are both also quite simple and effective. 
The additional switching electrode at the exit of the ionizer 
can be implemented with minor modifications and uses the same 
circuitry as the ionization switching. This approach 
overcomes the limitation imposed by the finite transit time 
down the ionizer. Suitable software can be written to 
measure and display average cluster size and distribution 
width, and this is discussed in additional detail below. The 
cluster size and distribution information is most valuable 
during beam set up and tuning prior to workpiece processing, 
but suitable circuitry, it is also possible to measure both 
the current and the cluster size during the overscan time of 
the mechanical scanning. 

Above it has been explained that gating off the ionized 
beam abruptly gives a signal in the Faraday cup that allows 
determination of the average mass of the clusters in the 
beam. The following discussion shows how to extract the full 
cluster number distribution from the signal. The cluster 
number distribution can be readily converted to the cluster 
mass distribution by multiplying by the mass of the molecules 
forming the clusters. Although the discussion refers to the 
measurement and processing of the beam shut off current fall 
signals when the beam is gated off, it is readily appreciated 
that a fully analogous situation occurs when ' the beam is 
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abruptly restored (gated on) and measurement and processing 
of the beam restoration current rise can equally be used as a 
means of implementing the invention. 
62 As previously explained, the switch shutting off the 

beam must be very fast compared to the ion flight times over 
the distance, L. Immediately after gating, the beam then 
consists of a column of ionized clusters of varying masses 
and velocities traveling from the gate to the Faraday cup. 
The beam current, I, received by the Faraday cup will begin 
to shut off when the smallest and, therefore, the fastest 
clusters travel the distance from the gating point to the 
Faraday cup. Current will be shut off completely when the 
largest (slowest) clusters travel this same distance. A 
number of techniques for accomplishing the switching were 
discussed previously. For simplicity, it is here assumed 
that the clusters are gated at the exit of the high voltage 
electrodes, so that the critical path length is the distance, 
L, between the gate at the exit of the high voltage 
electrodes and the Faraday cup. For a Gaussian number 
distribution the beam current in the Faraday is 
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Eqn. 5 



64 where N is the average number of molecules per cluster, a is 
the standard deviation, and 
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Here v(n) is the velocity of a cluster with n molecules, 
namely, 
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67 v(n)= 2q Vacc Eqn. 7 
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68 Note that a particular size cluster, n, does not affect 

current until 



69 VW t Eqn. 8 



70 At a given time, t, the cluster size that affects 

current is 
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72 where nip represents the atomic mass unit, ma represents the 
atomic mass number of the type of the individual atoms or 
molecules from which the cluster are formed, and n represents 
the number of molecules (atoms in the case of monatomic 
gases) in the cluster. 

73 Let us assume that there are so many different size 
clusters that we approximate the result for the current by a 
continuous rather than discrete variable. Then the current 
becomes . 
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Taking the negative derivative of this expression yields 
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79 Note that, nowhere in this derivation was assumption 

made that the distribution is Gaussian or any other form. 
Thus, more generally, the normalized probability distribution 
is 
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81 Note that the distribution must be plotted against n, as 
given in Eqn. 9. 

82 The foregoing derivation has assumed that the clusters 
all bear the same charge, q, a multiple of the of the unit 
charge, e. However, in general the clusters have a 
distribution in both size, n, and in charge, q. At the 
present, there is no easy separation of these distributions. 
Consequently in this case the time-of -flight measurement 
measures the joint distribution. Note that the velocity of 
the particles is dependent on n/q (Eqn. 8) and that the 
solution for n (Eqn. 9) may just as readily be expressed as a 
function of n/q. This is a result of the electrostatic 
nature of the acceleration system. In any event, all of the 
foregoing derivations follow in exactly the same manner as 
above except that the probability distribution function 
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becomes likewise f(n/q) rather than merely f(n). While this 
generalization somewhat reduces the utility of the 
measurement, the time-of -flight derived f(n/q) distribution 
is very useful in determining consistency of operation in 
commercial production equipment and for facilitating the set- 
up of equipment operating parameters to achieve consistent 
processing. Moreover, it is often possible to infer the size 
distribution separately by operating the cluster ionizer at 
low levels so that it is unlikely that the clusters acquire 
more than a single charge each - in such case the assumption 
that q is the same for all cluster ions is justified. 

A computer simulation has been written to perform the 
extraction of cluster size as outlined above by simulating 
the abrupt gating off of a cluster ion beam. The computer 
generates the current shut off signal based on the discrete 
sum (Eqn. 5), using either the Gaussian probability 
distribution above or a Poisson distribution of the form 



where m is the order of the distribution. 

The simulation is a discrete approximation that uses 50 
equally spaced points between the 100% and 0% points of the 
current shut off signal in exactly the same way that a 
periodically sampling data acquisition system does. The 
normalized derivative defined by Eqn. 13 is then calculated. 
When plotted against cluster size in the beam, the derived 
distribution is compared with the initial distribution. Some 
results of the simulation are shown in Figures 12 and 13 and 
in Figures 15 and 16 and discussed in below. 

Figure 7, as previously described, shows the current 
signal for a 20 keV beam with a cluster size distribution 
that is normally distributed and which has an average size of 
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1500 argon molecules per cluster and a standard deviation of 
400 molecules. 

88 Figure 12 shows the normalized negative derivative of 
the current signal from Figure 7 divided by time. The 
variations (wiggles in the curve) are the result of the; 
numerical approximation derivative in the simulation. 

89 Figure 13 shows this same distribution of Figure 12 
plotted against cluster size rather than time. The derived 
distribution is compared to the initial Gaussian 
distribution. Note that the agreement is excellent despite 
the numerical approximation of the derivative. 

q 90 Figure 14 shows the current shut off for a cluster size 

distribution having Poisson statistics of order 2 . Note that 
O in this case the average cluster size does not occur at the 

'fl 50% point. For the Poisson distributions, the average occurs 

(4 below the 50% point. As the order is increased then the 

average approaches the 50% point. 

91 Figure 15 shows the normalized time derivative of the 
signal of Figure 14 divided by time. It is recognized that 
smoothing techniques could be employed to improve the 
faithfulness of the derivative. For the Poisson 
distribution, the peak of the number distribution occurs at a 
number well below the average, so that more clusters have a 
lower time-of -flight than for the Gaussian distribution. 
This is offset by the longer tail at high masses, thus giving 
the same average . 

92 Figure 16 shows the normalized time derivative of the 
signal of Figure 14 plotted against the corresponding cluster 
size. Again, note the excellent fit, despite the rather 
dramatic change in the cluster distribution (Poisson versus 
Gaussian) . 

93 The process described here allows a full extraction of 
the cluster size distribution in the beam. Note that the 
process involves taking a derivative of the signal and 
dividing by time. The simulation illustrated above mimics 
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very well the deconvolution process that would be followed 
assuming the signal is taken as periodic samples with a 
conventional data acquisition system. Note that in 
comparison to the more usual pulsed time of flight, the 
signal levels here are high level and easily read with low 
noise. The simulation can easily be modified to take into 
account a finite switching time. This is expected to give a 
slight offset in the derived distribution but for the cluster 
sizes expected, the offset should be negligible. While the 
emphasis here has been on cluster statistics, any molecular 
ions present in the beam can also be detected. In a beamline 
including a molecular ion filter, none are ordinarily present 
and the possible failure of the molecular ion filter is 
detectable . 

Figures 17, 18, and 19 show results of measurements and 
analysis made using the invention in a gas cluster ion beam 
processing system as shown in Figure 4, where the beam gating 
device 212 is of the deflection switching type as previously 
detailed in Figure 11. A beam filter 214 was present in the 
gas cluster ion beam processing system, but was purposely 
disabled so that molecular ions and low mass ions would be 
present in the GCIB. An argon gas cluster ion beam 
accelerated by a 3kV acceleration voltage, V acc , was used for 
this example. The cluster ion beam current, I, was on the 
order of 10 microamps and was collected in a Faraday cup. 
The flight distance, L, was approximately 34 cm. The 
measurement system 226 collected periodic digitized samples 
of the collected beam current, I, and processed according to 
the principles of the invention discussed hereinabove. 

Figure 17 shows a graph of the measured gas cluster ion 
beam current, I, shut off characteristics when the beam was 
gated off by the deflector. A prompt current decrease of 
about 10% of the initial value of the current is noted near 
to time zero due to the short time-of- flight of the molecular 
ions and other very light ions. By 200 microseconds all the 
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beam current received in the Faraday cup is completely cut 
off. 

9 6 Figure 18 shows a graph of the time derivative of the 

measured current shut off signal shown in Figure 17 . 

97 Figure 19 shows a graph of the f (n/q) cluster 
probability distribution function corresponding to the 
current shut off characteristic measurement in Figure 17 and 
obtained according to the invention, plotted against n/q. 
Note the peak near n/q = zero, representing the molecular 
ions and smallest clusters and the later peak near n/q = 1300 
representing the peak of the main cluster distribution. This 
graphical information is produced by measurement system 22 6 
and displayed by display system 240 of the invention. 

98 Although the invention has been described with respect 
to various embodiments, it should be realized this invention 
is also capable of a wide variety of further and other 
embodiments within the spirit and scope of the appended 
claims . 

99 What is claimed is: 



